The infrared spectra of ABC-stacking tri- and tetra-layer graphenes 
studied by first-principles calculations 



O 

(N 

> 

O 

O 



c 

o 
o 



> 
in 
in 

ON 

d 



X 



Yuehua Xu 1 and San-Huang Ke'Q 

'Department of Physics, Tongji University, 1239 Siping Road, Shanghai 200092, P. R. of China 

The infrared absorption spectra of ABC-stacking tri- and tetra-layer graphenes are studied using the density 
functional theory. It is found that they exhibit very different characteristic peaks compared with those of AB- 
stacking ones, caused by the different stacking sequence and interlayer coupling. The anisotropy of the spectra 
with respect to the direction of the light electric field is significant. The spectra are more sensitive to the stacking 
number when the electric field is perpendicular to the graphene plane due to the interlayer polarization. The high 
sensitivities make it possible to identify the stacking sequence and stacking number of samples by comparing 
theory and experiment. 



I. INTRODUCTION 

Single-layer graphene (SG) is a two-dimensional (2D) flat 
monolayer composed of carbon atoms arranged in a honey- 
comb network, which is a basic building block for all other 
graphite materials. In the past, it was believed to be unstable 
in free-standing state [1] and so described as an 'academic' 
material. However, SG was unexpectedly made experimen- 
tally several years ago [2]. The unique electronic properties of 
SG is mainly due to its very peculiar band structure, with the it 
and 7r* bands showing linear dispersion around the Fermi level 
(E-p) where they touch with each other at a single point K in 
the Brillouin zone (BZ) [3, 4]. Because of its great potential 
for applications in nano-science and technology, SG has been 
thoroughly studied both theoretically and experimentally d5|— 

m. 

Recently, a further progress in this field was the successful 
fabrication and use of few-layer graphenes (FGs) which are 
stacking of a few graphene layers Ill| - [l3ll . From the view- 
point of application, FGs can be even more useful than SG 
since it offers further control of electronic states by interlayer 
interactions lfl4T - [l8ll , or by applying an electric field perpen- 
dicular to the molecular plane to open a band gap which is 
critical for applications in nanoelectronics H19L 12011 . Because 
of the interlayer interaction, theoretically, there can be dra- 
matic changes in the electronic properties of FGs compared 
with those of SG, depending on the stacking sequence and 
number (TV). Considering an arbitrary arrangement for adja- 
cent graphene layers, there will be 2 N ~ 2 possible low-energy 
stacking sequences for FGs of N layers. It is therefore de- 
sirable and important to understand the electronic properties 
of the different FGs and to be able to identify accurately their 
stacking sequence and number. Experimentally, infrared (IR) 
absorption spectra can be used to obtain the detailed informa- 
tion about the low-energy electronic excitations in FGs. If one 
has the knowledge about the influence of different stacking 
sequences and numbers on the IR absorption spectra, one can 
identify them accurately in terms of the experimental IR spec- 
tra. This knowledge can be obtained from theoretical calcula- 
tions, especially from first-principles calculations which have 
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the advantage of being empirical parameter free and therefore 
having the predictive power. 

Motivated by the fact that natural graphite adopts the AB- 
stacking (Bernal) sequence I2II1 which gives the lowest total 
energy, many theoretical and experimental studies in the past 
few years were focusing on the electronic and optical prop- 
erties of AB-stackingFGs I22H27I1 . especially their infrared 
absorption spectra Il28l - [3011 . For example, it was found that 
the interlayer coupling in an AB-stacking bilayer graphene 
(BG) leads to remarkable changes in its low-energy disper- 
sions as well as its IR spectra compared with those of SG 
lfl2[|3Tl - [34ll . On the other hand, very recently, ABC-stacking 
(rhombohedral) FGs were also found experimentally. For ex- 
ample, Norimatsu et al observed selective formation of ABC- 
stacking graphene layers JHI, and more recently, Mak et al re- 
ported unambiguous experimental evidence for the existence 
of stable tetra-layer graphenes in both AB- and ABC-stacking 
sequences 113611 . The existence of these stable polytypes of FG 
provides a new possible way to tailor the electronic structure 
of FGs materials for both fundamental studies and application 
interests [37-40]. However, because of the lack of practical 
samples and related experimental observation in the past, the- 
oretical studies on the band structure and optical properties, 
especially the IR absorption spectra, of the ABC-stacking se- 
quence are still very limited 1141 144511. except for some simpli- 
fied tight-binding (TB) modeling 1461 l47tl . 

In this paper, we study the electronic structures and 
IR absorption spectra of ABC-stacking tri- and tetra-layer 
graphenes (for simplicity, we call them ABC-3 and ABC-4, 
respectively) using first-principles density functional theory 
(DFT) calculations. The low-energy band dispersion and the 
formation of the characteristic peaks in their IR absorption 
spectra are analyzed in comparison with those of the AB- 
stacking ones (called AB-3 and AB-4, respectively), showing 
significant effects from the stacking sequence. The anisotropy 
of the spectra with respect to the direction of the light electric 
field is also investigated and is found to be remarkable. Our 
theoretical results for AB-4 and ABC-4 are in good agreement 
with an very recent experimental report IB6I1 . Compared to TB 
model calculations, the present calculation has the advantage 
of better transferability and predictive power since it is free of 
any empirical parameter. We show that, together with reliable 
experimental data, the theoretical calculation of IR spectra can 
provide a useful tool for identifying the stacking sequence and 
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even the stacking number of samples. The rest part of our pa- 
per is arranged as follows: In the next section (Section II) we 
give briefly the theory and computational details; We discuss 
the results in Section III, followed by a summary in Section 
IV . 



II. THEORY AND COMPUTATIONAL DETAILS 

A. IR optical absorption 

The ground-state electronic structure of a FG can be calcu- 
lated using density functional theory. After the band structure 
E(k, j) and corresponding eigenfunctions are obtained, its IR 
absorption properties can be studied by looking at the imagi- 
nary part of the frequency-dependent dielectric function 114811 : 



e 2 r (uj) = — - — lim — } 2w k 

x 5[£(k,ji)-£;(k,j2)-w], (1) 

where ui and q are the energy and momentum of photon, 
respectively; ft denotes the volume of the unit cell; e is the 
charge of electron and % is the weight of k-point k for the 
k-sampling; Indices ji= 1 ,2, N and j% = -1, -2, — N 
denote the conduction and valence subbands, respectively, 
counted from the Fermi level; Ujjc is the cell periodic part 
of the wavefunction of band j at k-point k . Under an electric 
field parallel to the graphene plane ( E || x or E || y) or per- 
pendicular to the graphene plane (E || z ) which is studied in 
our work, the electrons could be excited from the occupied va- 
lence 7r bands to the unoccupied conduction ir* bands. At zero 
temperature, only inter- n -band excitations can occur with the 
excitation energy cj cx = E(k, ji) — E(k, j'2). Because pho- 
ton's momentum is almost zero, the optical selection rules is 
Ak = 0. 



B. Computational details 

For band structure calculation, we use the density func- 
tional theory implemented in the plane-wave pseudopotential 
formalism [49] with the local density approximation (LDA) in 
the Ceperley-Alder version B50I1 for the electron exchange and 
correlation. The interaction between the ions and electrons is 
described by the highly accurate full-potential projected aug- 
mented wave (PAW) method 115 ll 15211 which can give a more 
accurate and reliable result than the ultrasoft pseudopotential. 
In our calculations, the 2s and 2p orbitals of the carbon atoms 
are treated as valence orbitals, and a large plane-wave cut- 
off of 500 eV is used throughout. A uniform grid larger than 
200 x 200 x 1 is used in the irreducible BZ for the k-sampling. 

A supercell geometry is constructed so that the tri- and 
tetra-layer graphenes are aligned in a hexagonal supercell with 
the closest distance between the adjacent FGs being at least 



10A along the stacking direction (z direction). The interlayer 
spacing is set initially to be the same as in graphite. The struc- 
tures of ABC-3 and AB-3 are shown in Figs. 1(a) and 1(b), re- 
spectively, where there exist two C atoms, denoted by A4 and 
Bi , in a primitive cell of each graphene layer i. The two poly- 
types can be obtained by displacing each adjacent graphene 
layer continuously in one direction for the ABC stacking, and 
alternatively in opposite direction for the AB stacking. We use 
a conjugated gradient method to optimize the atom's positions 
as well as the size of the supercell until the forces acting on all 
the atoms are less than 0.005 eV /A. With the obtained equi- 
librium structure, we then use a more dense set of k-points to 
calculate the electronic structure and the moment matrix ele- 
ments. Finally we use the Eq. (1) to study the IR absorption 
properties. 



III. RESULTS AND DISCUSSION 
A. Atomic structure and energy band structure 

The equilibrium C-C bond length a given by the DFT-LDA 
structural optimization calculation is 1 .41 A for both the ABC- 
and AB-stacking FGs, which is almost the same as that in bulk 
graphite. However, the interlayer spacing b is slightly differ- 
ent for the different stacking sequences and stacking numbers: 
3.31 and 3.32A for ABC-3 and ABC-4, respectively, and 3.32 
and 3.34A for AB-3 and AB-4, respectively. Our result of 
3.31 A for ABC-3 is in good agreement with an available ex- 
perimental value 3.335A [45], both are a bit smaller than the 
value 3.37A in ABC-stacking graphite JH]. 

The low-energy band dispersions of the four systems are 
given in Figs. 2 (a) - (d). Let's first look at the AB-stacking 
case (Figs. 2(a) and 2(c)). As is evident, the low-energy band 
structure consists of separate single-layer-like and bilayer-like 
bands [44] but with some degree of distortion due to the ad- 
ditional interlayer couplings. Meanwhile, the interlayer cou- 
plings reduce the symmetry between the valence and conduc- 
tion bands with respect to the Fermi level, and also cause a 
small overlap between them near Ep. For example, for AB- 
3, our calculation predicts different energy splittings of -483 
and +552 meV at the K point for the low-lying and up-lying 
bands, respectively. For AB-4, the above band splittings at 
the K point given by our calculation are -551meV, -21 lmeV, 
222meV, 597meV for the four bands (see Fig. 2(c)), respec- 
tively, while they are ±222 and ±596 meV from a TB model 
calculation [36]. We note that our band structure is asymmet- 
rical with respect to the Fermi level while the result of the TB 
model lf36Tl is symmetrical. Based upon the previous experi- 
mental IR absorption spectra of doped or gated AB-stacking 
FGs HHH^I!!!], one can deduce that the valence and conduc- 
tion bands are asymmetrical with respect to E-p. Our result 
for the AB stacking is consistent with the experimental data 
while the result given by the simplified TB model is not, in- 
dicating that our DFT-LDA band structure calculation is more 
reasonable than the TB model calculation. 

For ABC-3 and ABC-4 which is the focus of our study, one 
can see in Figs. 2(b) and 2(d) that they have strikingly differ- 



( a ) AB C -s tacking tri-layer graphene (b ) AB -s tacking tri-layer graphene 



FIG. 1: Structures of tri-layer graphenes: (a) ABC-stacking, (b) AB-stacking. The C-C bond length is denoted by a and the interlayer spacing 
is denoted by b. 
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FIG. 2: The low-energy band dispersions of the four systems studied: (a) AB-3, (b) ABC-3, (c) AB-4, and (d) ABC-4, as used in the text. 



ent 2D band structures from those of the AB-stacking ones, 
which reflect the underlying difference in crystallographic 
symmetry and interlayer coupling. First, unlike in the AB- 
stacking case, the band dispersions of the ABC-stacking FGs 
don't resemble those of single-layer and bilayer graphenes. 
The different interlayer interaction causes two intersections at 
the K point with the energy of +360meV and -370meV, re- 
spectively, which are absent in AB-3 and AB-4. Again, here 
the band structure is slightly asymmetrical with respect to the 
Fermi level while it is symmetrical in some TB model calcula- 
tions (±360meV) ll36l[37ll . Second, the lower crystallographic 



symmetry of the ABC-stacking shifts the extrema of the low- 
lying and up-lying energy bands away from the K point. 
In ABC-3 case, the four shifted extrema of these "wizard- 
hat" bands have energies of about 335meV away from Ep, 
as shown in Fig, 2(b), while in ABC-4 case they are about 
290meV away from Ep. Third, the interlayer interactions pro- 
duce one pair of localized flat bands near Ep. For ABC-3, the 
two flat bands have a tiny energy gap on the K — T line while 
for ABC-4 they have a small overlap. Physically, these flat 
bands are formed by the localized electronic states from the 
two outermost graphene layers. 
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FIG. 3: Calculated IR absorption spectra of (a) ABC-3 and AB-3 under E || x, (b) ABC-4 and AB-4 under E || x, (c) ABC-3 and AB-3 under 
E || z, and (d) ABC-4 and AB-4 under E || z. The solid line is for ABC stacking and the dotted line is for AB stacking. The corresponding 
joint densities of states (JDOS) are shown in (e) and (f). 



B. IR absorption spectra 

We calculate the IR absorption spectra of the AB- and 
ABC-stacking tri- and tetra-layer FGs for two different direc- 
tions of the light electric field (E || x and E || z) and plot 
the results in Fig. 3 (a) - (d) together with the correspond- 
ing joint densities of states (JDOS) in Fig. 3(e) and (f). Let 
us start with the discussion about the different behavior be- 
tween the ABC stacking (the solid lines in Fig. 3) and the 
AB stacking (the dotted lines). We first consider the case of 
E || x. As one can see in Fig. 3(a), for AB-3 there exists 
one characteristic peak at 555 meV, which is mainly caused 
by the transitions between the valence bands j 2 = -3 (-1) and 
conduction bands j\ =1 (3) with their average energy differ- 
ence being about 555 meV around the K point, denoted by 
the red arrows in Fig. 4(a). For ABC-3, on the other hand, its 
characteristic peak is located at 355 meV. This peak is due to 
the transitions between valence band j 2 = -2 (around the ex- 
tremum) and conduction band ji=l as well as j 2 =-\ and j\ = 
2 (around the extremum) with their average energy difference 
being 355meV, as denoted by the red arrows in Fig. 4(b). The 
corresponding peaks in the JDOS can be found in Fig. 3(e). 
For ABC-3 the peak at 355 meV in the JDOS is due to ID-like 
van Hove singularity (vHS) which was analyzed earlier by a 
TB model 56 and ascribed to the fact that the "wizard-hat" 
bands have their extrema away from the K point. In addition, 
for ABC-3, there also exists another weak peak at 750 meV, 
which is mainly caused by transitions between the two inter- 



sections around the K point, i.e., between the valence bands 
j 2 = -2, -3 and conduction bands j\ = 2, 3. Obviously, the 
characteristic peak at 355meV of ABC-3 has a red-shift of 
200meV compared with that at 555 meV of AB-3 while the 
peak at 750meV of ABC-3 is absent in AB-3. For AB-4, as 
one can see in Fig. 3(b), there are two characteristic peaks, 
one at 225 meV and the other at 620 meV. The peak at 225 
meV is mainly caused by the transitions between the valence 
bands j 2 = -3 (-3) and conduction bands ji =1 (2) with the 
interlayer-coupling induced splitting energy between the pairs 
being about 225 meV, which are denoted by the red arrows in 
Fig. 4(c). The peak at 620 meV is due to the transitions be- 
tween the valence bands j 2 = -2, -1, -4 and conductance bands 
ji = 4, 4, 1, respectively. For ABC-4, on the other hand, there 
is one major peak at 309 meV which is mainly caused by the 
transitions between the valence band j 2 = -2 (around the ex- 
tremum) and conduction band ji = 1 as well as j 2 = -1 and ji 
= 2 (around the extremum) with the average energy difference 
309 meV, denoted by the red arrows in Fig. 4(d). The cor- 
responding peak in the JDOS (ID-like vHS) can be observed 
at around 309 meV in Fig. 3(f). Besides the major charac- 
teristic peak there also exist two weaker structures in the IR 
spectra of ABC-4. One is a "shoulder structure" at 356 meV, 
which is mainly caused by the transitions between the valence 
bands j 2 = -3 (-1) and conduction bands j\ = 1 (3) around 
the K points. The other is a weaker peak at 705 meV which 
is induced by the transition between the intersections of the 
valence bands j 2 = -2 (-3) and conduction bands j\ = 3 (2) 
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FIG. 4: Schematic drawing of the band structures of (a) AB-3, (b) ABC-3, (c) AB-4, and (d) ABC-4, which are in the same order as in Fig. 2. 
Red and black arrows denote different permitted optical transitions. Conduction bands are denoted by index jj and valence bands are denoted 
by index fa ■ 



around the K points. Overall, the three structures in the IR 
spectrum of ABC-4 are blue-shifted with respect to the two 
structures of AB-4. Specifically, the major peak at 309 meV 
has a blue-shift of 84 meV with respect to the peak at 225 
meV of AB-4. The weak peak at 705 meV of ABC-4 has a 
blue-shift of 85 meV with respect to the peak at 620 meV of 
AB-4. While the "shoulder structure" at 356 meV is absent in 
AB-4. 

Next, we consider the case of E || z (see Figs. 3(c) and 
(d)). For AB-3, a jump lies at 538 meV, which is caused by 
transitions between valence bands fa = -2 (-3) and conduction 
bands fa = 3 (2) near the K point with their average energy 
difference being about 538 meV, denoted by the black arrows 
in Fig. 4(a). The corresponding jump structure in the JDOS 
can be seen in Fig. 3(e). In contrast, for ABC-3, a charac- 
teristic peak appears at 690 meV. The transitions behind are 
between the valence band fa = -2 (around the extremum) and 
conduction band fa = 2 (around the extremum) as well as fa 
= -3 and fa = 3, denoted by the black arrows in Fig. 4(b). 
This peak corresponds to the ID-like vHS at 690 meV in the 
JDOS (see Fig. 3(e)), which is quite different from the 2D 
systems' characteristic jump structure (or call step singular- 
ity) at 538 meV of AB-3. For AB-4, as can be seen in Fig. 
3 (d), there is a weak peak at 230 meV and a major peak at 
455 meV which is induced by the transition between fa = - 
3 and fa = 3, denoted by the black arrows in Fig. 4(c). In 
Fig. 3(f) one can see the corresponding features in the JDOS 
for these two peaks. For ABC-4, on the other hand, a major 



peak is at 595 meV, as shown in Fig. 3(d), which is due to the 
transitions between the extrema of valence band fa = -2 and 
conduction band fa = 2, denoted by the black arrows in Fig. 
4(d). The corresponding peak in the JDOS (ID-like vHS) is 
shown in Fig. 3(f). Another weaker peak at 760meV (Fig. 
3(d)) is caused by the transition between fa = -3 and fa = 3 
around the K point. Overall, the peaks in the IR spectrum of 
ABC-4 are blue-shifted with respect to those of AB-4. 

The different IR absorption behavior between the ABC- 
and AB-stacking FGs can be understood by considering the 
difference in their electronic structure. As mentioned previ- 
ously, the ABC stacking has a lower crystallographic sym- 
metry compared with the AB stacking. As a result, the ex- 
trema of its low- and up-lying energy bands are shifted away 
from the K point, which gives rise to the ID-like vHS in their 
JDOS and prominent optical absorption peaks. These promi- 
nent peaks are in strong contrast with the much weaker inten- 
sity structures in the IR spectra of the AB stacking, and also 
located at very different frequencies . 

So far, almost all experiments ||32[ 1361 Bil l55ll and theoret- 
ical TB model calculations lfl5i l46l l47ll considered only the 
case of E || x in investigating the IR spectra of FGs, while the 
case of E || z was little concerned [34] although the setup of 
E || z is feasible in experiment and can provide more infor- 
mation, as already implied in our previous discussion. Below 
we discuss in more details the effect of different directions 
of the light electric filed on the IR absorption spectra. It is 
evident from Fig. 3 that the IR absorption spectra show a 
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FIG. 5: IR optical conductivity spectra of ABC-4 (solid line) and 
AB-4 (dashed line) under E || x. The inset shows the corresponding 
experimental results reported by Mak et al in Ref.[36]. The positions 
of the peaks are denoted. 

strong anisotropy between E || x and E || z. For ABC-3, 
the position of its major characteristic peak (at 690 meV) un- 
der E || z has a large blue-shift of 335 meV compared with 
the case of E || x. In addition, there is a weak peak at 750 
meV under E || x, which is absent under E || z. For AB- 
3 the peak at 555 meV under E || x is replaced by a jump, 
instead of a peak, around 538 meV, as shown in Fig. 3(c). 
Similar to the tri-layer systems, the tetra-layer systems also 
show a strong anisotropy effect between E || x and E || z, 
particularly ABC-4: Its characteristic peaks are significantly 
blue-shifted from E || x to E || z, by up to 286 meV. One 
thing to note is that, compared with the AB stacking, the ABC 
stacking has a much stronger anisotropy effect causing large 
blue-shifts for its major peaks from E || x to E || z (see Fig. 
3 (a) vs. (c), and (b) vs. (d)). Physically, this can be under- 
stood also by considering its lower crystallographic symmetry 
which results in a very different band structure including the 
"wizard-hat" bands. This band structure provides two sets of 
transitions with significantly different transition energies for 
E || x and E || z, respectively (see Fig. 4, red arrows vs. 
black arrows). 

Finally, we would like to discuss how the number of stack- 
ing layers affects the IR spectra. As can be seen in Fig. 3, the 
thickness effect ((a) vs. (b) and (c) vs. (d)) is not so signifi- 
cant compared to the effect from stacking sequence. For the 
ABC stacking, all the characteristic peaks have a red-shift as 
the stacking number increases. For the AB stacking, however, 
the trend is not so clear. Under E || x the maximum shift 
in peak position is about 45meV between ABC-3 and ABC- 
4. This shift is enhanced to more than 90meV under E || z. 
Physically, this enhanced effect can be understood by consid- 
ering the interaction between the light electric field and the 
graphene sheets. In the former case the field will induce an 
intra-layer polarization which is less sensitive to the thickness 
while in the latter the field will induce an inter-layer polariza- 
tion which is more sensitive to the thickness. 

From our discussion, one can see that different stacking se- 
quences of FGs have very significant effects on their IR ab- 



sorption spectra, leading to peaks at different frequencies and 
with different intensities. Furthermore, our calculation shows 
that the IR spectra are also sensitive to the direction of the 
light electric field, especially for the ABC stacking sequence. 
In the case of E || z the thickness of an ABC-stacking sample 
also plays an important role. The high sensitivities of the IR 
absorption spectra to the different factors provide important 
information for identifying the stacking sequence and stack- 
ing number of an experimental sample by comparing the ex- 
perimental spectra with the calculated ones. Below, we com- 
pare our theoretical result with an experimental report lr36Tl 
which is, to the best of our knowledge, the only available one 
in literature for ABC-stacking FGs so far. In order to have 
a reasonable comparison, we plot the IR optical conductivity 
spectra <t{uj) = we 2 (w)/47r Jlj] in Fig. 5 for ABC-4 and AB- 
4 under E || x. One can see that our theoretical results are 
in good overall agreement with the experiment data (see the 
inset in Fig. 5) for both ABC-4 and AB-4. The positions of 
the peaks and their relative intensities as well as the overall 
shape of the two experimental spectra are all well reproduced. 
Quantitatively, however, there are still some discrepancies be- 
tween theory and experiment. For AB-4 the maximum differ- 
ence in the peak position is about 45 meV while for ABC-4 
it is up to 50 meV. Additionally, for ABC-4 our calculation 
predicts a weak peak at 705 meV, which was also found in the 
simplified TB's model calculation [36] but with a slight shift 
(670 meV). However, this peak was not observed in the exper- 
iment. How to understand these quantitative discrepancies is 
still an open problem. It may be due to some experimental en- 
vironment which has not been taken into account in our DFT 
calculation, such as doping, disorder, defect effects. The ex- 
perimentally missed peak at 705 meV, on the other hand, may 
be due to temperature broadening. To clarify the quantitative 
inconsistency between theory and experiment, further theoret- 
ical work and more experimental measurements are desirable. 



IV. SUMMARY 

We have studied the IR absorption spectra of ABC-stacking 
tri- and tetra-layer FGs using a first-principle method for two 
different directions of the light electric field (E || x and 
E || z), and compare them with those of AB-stacking ones. 
The findings are as follows. 1) The ABC-stacking sequence 
causes great different band structures, inducing characteristic 
peaks at different positions and with different intensities, no 
matter E || x or E || z. 2) The IR spectra show a strong 
anisotropy effect between E || x and E || z. 3) The IR 
spectra are much more sensitive to the stacking number un- 
der E || z because of the induced inter-layer polarization. 4) 
Our calculated IR optical conductivity spectra for the tetra- 
layer graphenes under E || x are well consistent with a recent 
experimental observation. The significant effects of the differ- 
ent factors on the IR spectra of FGs provide useful tools for 
identifying their stacking sequence and stacking number. 

This work was supported by Shanghai Pujiang Pro- 
gram under Grant 10PJ1410000 and the MOST 973 Project 
2011CB922200. 
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